Multifractals of Normalized First Passage Time in 
Sierpinski Gasket 



Kyungsik Kim and J. S. Choi 
Department of Physics, Pukyong National University, Pusan 608-737, 

Korea 

Y.S. Kongf 

School of Ocean Engineering, Pukyong National University, Pusan 

608-737, Korea 



The multifractal behavior of the normalized first passage time is inves- 
tigated on the two dimensional Sierpinski gasket with both absorbing 
and reflecting barriers. The normalized first passage time for Sinai 
model and the logistic model to arrive at the absorbing barrier after 
starting from an arbitrary site, especially obtained by the calcula- 
tion via the Monte Carlo simulation, is discussed numerically. The 
generalized dimension and the spectrum are also estimated from the 
distribution of the normalized first passage time, and compared with 
the results on the finitely square lattice. 



PACS numbers: 05.45. +b, 05.60+w, 

Key words : normalied first passage time, generalized dimension, 
logistc map, Sierpinski gasket, multifractal 

E-mail: kskim@dolphin.pknu.ac.kr, jschoi@dolphin.pknu.ac.kr 
f E-mail : y skong® dolphin . pknu . ac . kr 



1 



1 Introduction 



Recently a lot of interest has been concentrated on the behavior of 
the disorder system in a variety of contexts in condensed matter 
physicfli. In the transport process one of the important subjects 
for the regular and disorder systems is the random walk theory, 
and this theory has extensively been developed to the continuous- 
time random walk theory which is essentially described both by the 
transition probability dependent of the length between steps and 
by the distribution of the pausing times i-i. Until quite recent the 
transport phenomena for the motion of a particle have been largely 
extended to the reaction kinetics iH, and the strange kineticS. 

One of the well-known interesting issues in the theory of the 
random walk is the mean first passage time that is defined by the 
average time arriving at the absorbing barrier for first time after 
a particle initially starts from an arbitrary site. In particular, this 
problem has extensively been studied in the Sinai model which is 
discussed on the random barrier model such as the absorbing and 
reflecting barriers. In previous work the Sinai model with asymmet- 
ric transition probability also was studied for the mean and mean 
square displacements dependent anomalously on time i-i. More- 
over, the transport process for the mean first passage time has 
already been addressed by Noskowicz et a 10 who have obtained 
the upper and lower bounds using the recursion - relation proce- 
dure™. In one dimensional lattice with the periodic boundary 
condition, Kozak et al@ treated long ago with the first passage 
time of a particle for the case of the transition probability as a 
chaotic orbits of logistic map0, and obtained the statistical value 
of the mean first passage time with the anomalous exponent. 

On the other hand, since the multifractal characteristics is di- 
rectly related to analyze the statistical property of the normalized 
first passage time, the multifractal for the mean first passage time 
has been studied intensively in connection with the random walk 
problem in the one dimensional lattice with both absorbing and 
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reflecting barriers^". In general, the multifractal investigation for 
the chaotic and disorder systems have been applied analytically 
and numerically to the transport phenomena 11^1123 and the random 
fractal structure having the time-dependent random potential in 
the diffusive motiorHH. Murphy et a@ recently has argued on 
the multifractal behavior for the mean first passage time in one 
dimensional lattice having both absorbing and reflecting barriers. 
Therefore, we think that it is of great interest for the multifractal 
investigation to treat with the transport process described by the 
transition probability such as logistic map. 

The purpose of this paper is to investigate the multifractal be- 
havior of a particle executed on the random walk in the two di- 
mensional Sierpinski gasket. Specifically, we deal with the random 
walk for one particle having the transition probability given by the 
logistic map, and concentrate on the normalized first passage time 
for the random walk with symmetric and asymmetric transition 
probabilities on Sierpinski gasket existing both absorbing and re- 
flecting barriers. In section 2, we present the transition probability 
expressed in terms of the logistic map, and the convenient formulas 
of the normalized first passage time and the multifractal are intro- 
duced. In section 3, the distribution of normalized first passage 
time arriving at an absorbing barrier is considered in both Sinai 
model and the logistic model which has the transition probability 
given by logistic map. For simplicity, the generalized dimension 
and the spectrum from the statistical value of normalized first pas- 
sage time are also calculated numerically, and compared with the 
results reported earilierll!. Lastly, we give some conclusions. 



2 Normalized first passage time and multifractals 

First of all, we present the transition probability as a chaotic orbits 
of logistic map in this section. We also introduce the formulas of 
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the normalized first passage time, the generalized dimension, and 
the spectrum. 

Next, we shall consider the random walk of a particle in two- 
dimensional Sierpinski gasket in which the number of sites N n hav- 
ing stage n in the d-dimensional Sierpinski KaskeH is given by 
N n = (d + 1)(1 + (d + l) n )/2. It is assumed that a particle is 
started from an arbitrary site A on Sierpinski gasket. The re- 
flecting barriers are located on all the boundary, except that an 
absorbing barrier is at an ended site B under the right side, as 
shown in Fig.l. 

Now, the transition probability for the motion of a particle is 
introduced as 

x(n + 1) = Rx(n)[l — x(n)] (1) 



where R is the controll parameter with < R < 4. In Monte 
Carlo simulation one has to treat with the arbitrary values of R, 
where the numerical generated sequences appear to be chati<E2lo. 
By using the transition probability as Eq. p], one particle jumps 
to right site with pij + 7 or left site with q\j — 7 , and to up site 
with p 2 j + 7 or down site with q 2 j — 7 in four directions after a 
particle starts from a site j on two-dimensional Sierpinski gasket, 
where 7 is called the disorder parameter. It has really been showed 
that the normalized condition for transition probability is given by 
Pij + Qij + P2j + <l2j = 1- In particular, the transition probability 
is symmetric for the disorder parameter 7 = 0, and for < 7 < j, 
it also assumes that a particle executes on a biased random walk 
jumped to the direction of the absorbing barrier. 

Using the generating function technique the mean first passage 
timJIII < T > in one dimensional lattice is given by 

N-l I N-2 1 AT-l i 

<^>=E^ + E^ E n I- P) 

k=l yk k=0 fk i=k+l j=k+l t>3 



where one particle moves to right site with pj or left site with gr- 
after one step, when it starts from site j. However, the transi- 
tion probability has always the different value at each site when a 
particle jumps for ocurrence of the one step to a nearest-neighbor 
site from a given site, and the derivation for Eq. ^ can be found 
in details elsewheriJli3t3. Particularly, a possible extension of this 
paper is only to two dimension in order to discuss briefly on the 
mean first passage time and the multifractal in two dimensional 
Sierpinski gasket. The normalized first passage time T nt has the 
form 

rp t-n t-min /q\ 
J-nt — ~ — \0) 

where t n is the first passage time arrivng at the absorbing barrier 
after n steps, and t max and t m i n are the maximum and minimum 
values of the first passage time arrivng at the absorbing barrier, 
respectively. In the next section we will perform the conventional 
numerical simulations for the normalized first passage time in both 
Sinai model and the logistic model on two-dimensional Sierpinski 
gasket. 

To show the multifractal feature for the motion of a particle 
with the transition probability expressed in terms of the logistic 
function, we extend to the multifractal calculations for the normal- 
ized first passage time. If we divide the normalized first passage 
time into a set e of as e — > 0, then the generalized dimension in the 
multifractal structurSS is represented as 



lnJ2i riip q - 
**6 (q -l)lne 



D q = lira V_ ; ( 4 ) 



where Pi is the probability for the number of configuration of i 
particle arrived at the absorbing barrier and rij the number of con- 
figuration for i th particle. By introducing the above expression the 
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spectrum f q and a q is caculated from the relations 

f q = a q -(q- l)D q 



(5) 



and 

« 9 = ^[(?-l)4 (6) 
where Eq. [5] and Eq. ^| can be obtained by Legendre transformation. 
3 Calculation and results 

We mainly concentrate on the generalized dimension and the spec- 
trum for the normalized first passage time on two-dimensional Sier- 
pinski gasket lattice. We only restrict ourselves to the case of n = 7 
stages on Sierpinski gasket, even though the stage n can be ex- 
tended to large number in this paper. At first, three sites are at 
(1), (2), and (3) in the initial stage n — 0, and we assume that 
a particle starts from an initial point iq = (11111233) at n = 7 
stages on two-dimensional Sierpinski gasket lattice. The reflecting 
barriers is located on all the boundary surface, except that an ab- 
sorbing barrier is at B = (33333333) in Fig.l. Next, both the cases 
of Sinai model and the logistic model, for asymmetric as well as for 
symmetric transition probability, are considered as follows: One is 
the case of Sinai model in which it has the symmetric transition 
probability with = qij = p 2 j = q 2 j = \ and 7 = 0, and the 
two asymmetric transition probabilities with 7 = 0.05 and 7 = 0.1. 
The other is the case of the logistic model which is related to the 
transition probability given by the logistic map under the same 
condition of Sinai model. 

From now on, we estimate numerically the generalized dimen- 
sion and the spectrum after finding the normalized first passage 
time from Eq. |3] via Monte Carlo simulation. To discuss on these 
multifractal quantities for the normalized first passage time in two 
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dimensional Sierpinski gasket, we perform the numerical simula- 
tions to show our analytic arguments. For three values of the dis- 
order parameter (i.e., 7 = 0, 0.05, 0.1) in both Sinai model and the 
logistic model, our simulations are performed 3 x 10 6 particles and 
averaged over 10 4 configurations, and analyzed the normalized first 
passage time, the generalized dimension, and the spectrum. The 
result of these calculations is summarized in Table 1. 

Particularly, in the logistic model we are actually interested in 
the one case for the control parameter R = 3.9999. It can be easily 
found from the iteration of the logistic map that three symmet- 
ric cutoff values are 0.152210, 0.505753, and 0.855890 for 7 = 0, 
while three asymmetric cutoff values are 0.280024, 0.505753, and 
0.730526 for 7 = 0.1, respectively. Here the cutoff value is defined 
as the quantity used to determine the direction of random walker 
in the logistic model. 

As shown in Table 1, the fractal dimenion D Q ( i.e., the maxi- 
mum value of D q or f q ), the scaling exponent a q , and the general- 
ized dimension D q , ultimately based on the theoretical expressions 
Eq. [| - Eq. are calculated numerically in both Sinai model and 
the logistic model. In Sinai model Figs. 2 and 3 depict respectively 
the generalized dimension D q as a function of q and the spectrum 
f q as a function of a q for varying the disorder parameter values. 
It can be seen from Table 1 that the fractal dimension changes 
anomalously as the disorder parameter 7 increases in two models. 
In particular, it is also found from the result obtained in our sim- 
ulations that in Sinai model for 7 = 0.05 the value of the fractal 
dimension on Sierpinski gasket is nearly equal to that on square 
lattice. Futhermore, in the logistic model on the two dimensional 
Sierpinski gasket, the normalized first passage time is found to be 
infinite for 7 = 0, and the fractal dimension is expected to take the 
value near zero as 7 — > 0.5. 

In conclusion, we have here considered the normalized first pas- 
sage time for both Sinai model and the logistic model on two dimen- 
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sional Sierpinski gasket. We have clearly showed the multifractal 
characteristics from the distribution of the normalized first passage 
time, as summarized in Table 1. In future work, for the multifractal 
characteristics for the normalized first passage time, we intend to 
compare with the results found in both Sinai model and the logistic 
model if the transition probability with the modified log - normal 
function is introduced, and to attempt to investigate extensively in 
the similar lattice models. 
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Figures 



Fig.l. The trajectory of a particle in the two dimensional Sierpi- 
inski gasket at n = 7 stages. A particle is started from at 
A = (11111233), and the absorbing barrier represents the 
site B = (33333333). 

Fig. 2. The generalized dimension D q versus q for the normal- 
ized first passage time on the two dimensional Sierpinski 
gasket in Sinai model. 

Fig. 3. The spectrum f q versus a q for the normalized first pas- 
sage time on the two dimensional Sierpinski gasket in Sinai 
model. 



Table 



Table I. Values of the fractal dimension, the generalized dimen- 
sion, and the scaling exponent in Sinai model and the logis- 
tic model. These lattice models are on both Sierpinski 
gasket and the finite square latticH. 
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